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Abstract 

Let A = {ai, a 2 , ■ ■ ■ , a s } and {bi, b 2 , ■ ■ ■ , b t } be two sets of objects with 
s + r = n, the generalized assignment problem assigns each element cij G A 
to at least ctj and at most a'i elements in B, and each element bj G B to 
at least /3j and at most • elements in A for all 1 < % < s and 1 < j < t. 
In this paper, we present an 0(n 4 ) time and 0(n) space algorithm for this 
problem using the well known Hungarian algorithm. We also present an 
0(n 3 ) algorithm for a special case of the generalized assignment, called the 
limited- capacity assignment problem, where a^, (3j = 1 for all i, j. 
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1. Introduction 

Given two sets A and B, the assignment problem aims to optimally assign 
each object of one set to at least one object of the other set. The assignment 
problem has applications in various fields such as computational biology 1], 
pattern recognition 0|, computer vision j^j, music information retrieval 4], 
and computational music theory |5|. Let A and B be two sets with + = 
n, Eiter and Mannila 0] proposed an 0(n 3 ) algorithm for the assignment 
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problem between A and B by reducing it to the minimum-weight perfect 
matching problem in a bipartite graph. 

In this paper, we consider the generalized assignment between A and B, 
called GA problem, where each point ai G A is assigned to at least ati and at 
most a'i points in B, and each point bj G B is assigned to at least /3j and at 
most (3'j points in A, such that sum of the matching costs is minimized. We 
also present an 0(n 3 ) time algorithm for a special case of the GA problem, 
where each object must be matched to at least one object. Schrijver [7J solved 
the GA problem in strongly polynomial time. We present a new algorithm 
which computes a generalized assignment between A and B in 0(n A ) time 
using 0(n) space. In section we review the basic Hungarian algorithm and 
some preliminary definitions. In section [3J we present our new algorithms. 

2. Preliminaries 

Given an undirected bipartite graph G — (AU B, E), a matching in G is 
a subset of the edges MCE, such that each vertex v G A U B is incident 
to at most one edge of M. Let Weight (a, b) denote the weight of the edge 
(a, b), the weight of the matching M is the sum of the weights of all edges in 
M, hence 

Weight(M) = ^ Weight(e). 

eeM 

A max-weight matching M is a matching that for any other matching M', 
we have Weight(M') < Weight(M). 

A path with the edges alternating between M and E — M is called an 
alternating path. Each vertex v that is incident to an edge in M is called 
a matched vertex; otherwise it is a free vertex. An alternating path that 
its both endpoints are free is called an augmenting path. Note that if the 
M edges of an augmenting path is replaced with the E — M ones, its size 
increases by 1. 

A vertex labeling is a function I : V — > R that assigns a label to each 
vertex v G V. A vertex labeling that in which 1(a) + 1(b) > Weight(a, b) for 
all a G A and b G B is called a feasible labeling. The equality graph of a 
feasible labeling I is a graph G = (V, E{) such that Ei = {(a, b)\l(a) + 1(b) = 
Weight(a,b)}. The neighbors of a vertex u G V is defined as Ni(u) = 
{v\(v,u) G Ei}. Consider a set of the vertices S C V, the neighbors of S is 
N l (S) = [j ueS N l (u). 



Lemma 1. Consider a feasible labeling I of an undirected bipartite graph 
G = (A U B, E) andS CA with T = N^S) ^ B, let 

ai = min {/(cij) + l(bj) — Weighty, bj)}. 

If the labels of the vertices of G is updated such that: 

(l(v) — ol\ if v G S 
l(v) + ai ifveT , 
l(v) Otherwise 

then I' is also a feasible labeling such that E\ C E[. 

Proof. Note that / is a feasible labeling, so we have 
1(a) + 1(b) > Weight(a,b) for each edge (a,b) of E. 
After the update four cases arise: 

• a E S and b G T. In this case 

l'(a) + l'(b) = 1(a) -ai + 1(b) + a t = 1(a) + 1(b) > Weight(a, b). 

• a ^ S and b ^ T. We have 

I' (a) + l'(b) = 1(a) + 1(b) > Weight(a, b). 

• a ^ S and b G T. We see that 

I' (a) + l'(b) = 1(a) + 1(b) + ai > 1(a) + 1(b) > Weight(a, b). 

• a G S and b T. In this situation we have 

l'(a) + l'(b) = 1(a) - on + 1(b). 

Two cases arises: 

- 1(a) + 1(b) - Weight(a, b) = a,. So 

l'(a)+l'(b) = l(a)-ai+l(b) = l(a)-l(a)-l(b)+Weight(a,b)+l(b) = Weight(a,b). 
Hence, E t C Ey. 



— 1(a) + 1(b) — Weight(a,b) > a\. Obviously 

I' (a) + l'(b) = 1(a) -a t + 1(b) > Weight(a, b). 

□ 

Theorem 1. If I is a feasible labeling and M is a Perfect matching in E\, 
then M is a max-weight matching f^j. 

Proof. Suppose that M' is a perfect matching in G, since each vertex is 
incident to exactly one edge of M' we have: 

Weight(M')= Weight(a,b) < l ( v )- 

(a,b)eM' ve(AuB) 

(aub) K v ) * s an u PP er bound for each perfect matching. Now assume 
that M is a perfect matching in Ef. 

Weight(M) = Y,Kz)= E 

egM ve(AuB) 

It is obvious that M is an optimal matching. □ 
In the following, we briefly describe the basic Hungarian algorithm which 
computes a max-weight perfect matching in an undirected bipartite graph 
G = (A U B,E) with \A\ = \B\ = n. It is obvious that for computing the 
minimum cost many to many matching using the Hungarian algorithm we 
must weight each edge (a h b 3 ) by Weigh ) {a ^y 

In lines 2 and 3, we label all vertices of B with zero and each vertex 
dj G A with max 7 j =1 Weight(a i ,bj) to get an initial feasible labeling. Note 
that M can be empty. In each iteration of the while loop of lines 5 — 21, two 
free nodes a; and bj are matched, so it iterates 0(n) times. Using the array 
skack[l..n], we can run each iteration of this loop in 0(n 2 ) time. The repeat 
loop runs at most 0(n) times until finding a free node bj. In line 11, we can 
compute the value of a\ by: 



oci = min slac, 



Algorithm 1 The Basic Hungarian algorithm^, £?) 



1: Initial > Find an initial feasible labeling I and a matching M in E\ 

2: Let l(bj) = 0, for all 1 < j <t 

3: /(ctj) = max* =1 Weight(ai, bj) for all 1 < % < s 

4: M = 

5: while M is not perfect do 

6: Select a free vertex ai E A and set 5 = {a«}, T = 

7: for j ; «— 1, n do 

8: s/acA;[j] = Z(aj) + — Weight^, bj) 

9: repeat 

10: if N t (S) =T then 

11: = mirifj^T slack [j] 

12: Update(l) > Update the labels according to Lemma [1] 

13: for all bj i T do 

14: s/acA;[j] = s/acA;[j] — on 

15: Select ue Ni(S) —T 

16: if m is not free then > (u is matched to a vertex z, extend the 
alternating tree) 

17: S = SU{z},T = TU{u}. 

18: for j •<— 1, n do 

19: sZacfcfj] = min(Z(2:) + l(bj) — Weight(z, bj), slack[j]) 

20: until u is free 

21: Augment(M) 
return M 



in 0(n) time. After computing ai and updating the labels of the vertices, 
we must also update the values of the slacks. This can be done using: 

for all bj ^ T, slack[j] = slack[j] — 

In line 12, we update the feasible labeling I such that Ni(S) 7^ T. In line 17 
of Algorithm (TJ when a vertex is moved form S to S the values of skack[l..n] 
must be updated. This is done in 0(n) time. 0(n) vertices are moved from 
S to 5*, so it takes the total time of 0{n 2 ). 

The value of ct\ may be computed at most O (n) times in O (n) , so running 
each iteration takes at most 0(n 2 ) time. So, the time complexity of the basic 
Hungarian algorithm is 0(n 3 ). 

3. The generalized assignment algorithm 

In this section, we describe our new algorithm which is based on the 
well known Hungarian algorithm. Consider two sets A = {ai, a 2 , . . . , a s } 
and B = {b\, 62, • • • , h} with s + t = n. Let Da = {«i, 012, ■ ■ ■ > «s} and 
Db = {Pi, 02, ■ ■ ■ , ftt} denote the demand sets of A and B, respectively. 
Let Ca = {a'l, ot'2, ■ ■ ■ , o! s\ and Cb = fl'21 ■ ■ ■ 1 ft't} be the capacity 
sets of A and B, respectively. Without loss of generality, we assume that 

The input of our algorithm is the complete bipartite graph that is con- 
structed as follows. Consider the complete bipartite graph G = (X U Y, E) 
where X = A U A' and Y = B U B' (see Figure [1]). A complete connection 
between two sets is a connection that in which each element of one set is con- 
nected to all elements of the other set. We show each set of the vertices by 
a rectangle and the complete connection between them by a line connecting 
the two corresponding rectangles. 

Given A = {a\, a 2 , . . . , a s } and B = {bi, 6 2 , • • • , h}, there exists a com- 
plete connection between A and B such that the weight of (tZj, bj) is equal to 
the cost of matching the point to bj for all 1 < i < s and 1 < j < t. Let 
B' = {b[,b' 2 ,...,b' t } and A' = , a' 2 , . . . , a' s } , each point of A is connected 
to the all points of B' such that the weight of (a^, b'j) is equal to the weight 
of (cii,bj). There exists also a complete connection between the sets B and 
A' such that the weight of (ct^, bj) is equal to the weight of (aj, bj). 

3.1. The generalized assignment algorithm 

Theorem 2. Let A and B be two sets with \A\ + \B\ = n, a generalized 
assignment between A and B can be computed in 0(n 4 ) time. 
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Figure 1: Our constructed complete bipartite graph with ^i=i a 'i — ELi P'j- 



Proof. We apply our new algorithm, Algorithm [21 on our bipartite graph 
G. Let Cap(u) and Dem(u) denote the capacity and the demand of the vertex 
u; so for all i,j we have Dem(ai) = a*, Dem(bj) = (3j, Cap{ai) = a-, and 
Cap{b j )=P' j . 

In our algorithm, a vertex x is free to another vertex y when x is not 
matched with y in M and has at least one empty capacity. So, E A and 
a- G A' are called free vertices to a vertex b that are not matched with it in 
M, if Num(ai) < Dem(a,i) and Num(a' i ) < Cap(ai) — Dem(ai), respectively. 
Also the vertices bj and b'j are free to another vertex that is not incident in 
M to them, when 

Num{bj) < Dem(bj) and Num{b'j) < Cap(bj) — Dem(bj), respectively. 

We save the current number of the vertices that are matched to the ver- 
tices of A, B, A', and B' in the arrays A[l . . . s], B[l...t], A'[l ...s], and 
B'[l . . . t], respectively; for example A[i] shows the number of the nodes that 
are matched to a*. The initial values of the arrays is 0; when a new point is 
matched to their representing node their values are increased by 1. Assume 
that Num(u) returns the number of the vertices that are matched to u so 
far. So Num(ai) = A[i], Num^a'^ = A'[i], Num(bj) = B[j], and finally 
NuTn(bj) = B'[j}. Note that the procedures IsFree(u) and IsMatched(u) 
return True if u is free and u is matched, respectively. Therefor, for example 



in the augmenting path a, 6, c, d, the vertex a is free to 6, b is matched to c, 
and d is free to c. Now we change the basic Hungarian algorithm as follows. 



Algorithm 2 The generalized assignment algorithm 
1: Initialize > Find an initial feasible labeling / and a matching M in E\ 



2: Let l{bj), Ity) = 0, for all 1 < j < t 

3: l(a,i) = max t j =1 (max(W eighty, bj), Weight^, b'^j) for all 1 < i < 
s 

4: K a 'i) = m ax* =1 Weight^, bj) for all 1 <i < s 

5: Let M = 

6: while {w G B U B' , with IsFree(w)} ^ do 

7: Select it G A U A' with IsFree(u) 

8: Set S = {w},T = 

9: for j <- l,t do 

10: sZacfcfj] = mm(l(u) + l(bj) — Weight(u, bj)) 

11: slack[j + t} = min(/(w) + l(b'j) - Weight(u, b'j)) 

12: repeat 

13: if A^(S) = T then 

14: cty = mi^min^x s/ac/c[j], min 6 ^ T s/ac/c[j + t\) 

15: Update(l) 

16: Select y e iV,^) -T 

17: if IsMatched(y) then 

> (y is matched to some vertices z) 

18: for all {z\(z,y) G M and z <£ S} do 

19: 5U{z},T = TU{i/} 

20: for j <— l,t do 

21: sZacfcfj] = min(/(z) + Z(6j) — Weight(z, bj), slack[j]) 

22: sZacfcfj' + i] = mm(l(z) + l(b' j )-Weight(z,b' j ),slack[j + 
t}) 

23: until IsFree(y) 

24: Augment(M) 



We first label the vertices of our bipartite graph G using an initial feasible 
labeling in lines 2 — 4. In each iteration of our algorithm, \M\ increases by 



I. So, our algorithm has 0(n 2 ) iterations with 0(n 2 ) time and runs in 0(n 4 ) 
time. 

□ 

3.2. The limited- capacity assignment algorithm 

Now we present an 0(n 3 ) algorithm for the limited-capacity assignment 
problem, where each object must be assigned to at least one point of the 
other set and the capacity of each object is limited. 

Theorem 3. Let A and B be two sets with \A\ + \B\ = n, a limited- capacity 
assignment between A and B can be computed in 0(n 3 ) time. 

Proof. We use the bipartite complete graph that is constructed for the 
generalized assignment problem. We modify the GA algorithm as following. 



Algorithm 3 The limited-capacity assignment algorithm (Part I) 
1: Initialize > Find an initial feasible labeling / and a matching M in E\ 



2: Let l(bj), 1%) = 0, for all l<j <t 

3: l(a,i) = max t j =1 (max(Weight(ai, bj), Weighty, bj)) for all 1 < i < 
s 

4: K a 'i) = m ax* =1 Weight^, bj) for all 1 <i < s 

5: Let M = 

6: while {u G A, with IsFree(u)} ^ do 

7: Select u E A with IsFree(u) 

8: Set S = {u},T = 

9: for j <- l,t do 

10: slack[j] = min(/(w) + l(bj) — Weight(u, bj)) 

11: s/acA;[j + t] = min(Z(u) + l(b'j) - Weight(u, b'j)) 

12: repeat 

13: if jVj(S') = T then 

14: = min(min fej ^ T slack[j], min 6 ^ T slack[j + t\) 

15: Update(l) 

16: Select y G JV,(5') - T 

17: if IsMatched(y) then 

> (y is matched to some vertices z) 

18: for all {z\(z,y) G M and z S} do 

19: 5U{z},T = TU{!/} 

20: for j <- l,t do 

21: s/ac/c[j] = mm(l(z) + l(bj) — Weight(z, bj), slack[j}) 

22: slack[j + t] =mm(l(z) + l(b'j)-Weight(z,b'j),slack[j + 



t}) 



23: until IsFree(y) 
24: Augment(M) 



Algorithm 4 The limited-capacity assignment algorithm (Part II) 



25: while {u e B, with IsFree(u)} ^ do 

26: Select u E B with IsFree(u) 

27: Set S= {«},T = 

28: for i 4- 1, s do 

29: slack[i] = min(/(-u) + /(aj) — Weighty, u)) 

30: s/ac/c[i + s] = mm(l(u) + Z(aQ — Weight^, u)) 

31: repeat 

32: if N t (S) — T then 

33: «z = min(min a .£ T sZacA;[i], min a ^ T s/ac/cfi + s]) 

34: Update(l) 

35: Select y e Ni(S) - T 

36: if IsMatched(y) then 

> (y is matched to some vertices z) 
37: for all {^|(^,y) G M and z ^ 5} do 

38: SU{z},T = TU{y} 

39: for i ■<— 1, s do 

40: s/ac/c[i] = min(/(z) + /(aj) — Weight(a,i, z), slack[i]) 

41: s/ac/c[i + s] — mm(l(z) + l(a f i ) —Weight(a' i ,u),slack[i + 

s}) 

42: until IsFree(y) 

43: Augment (M) 



□ 

4. Conclusion 

In this paper, we presented an 0(n 4 ) time and 0(n) space algorithm for 
computing a generalized assignment between A and B with total cardinality 
n. In fact, we modified the basic Hungarian algorithm to get a new algorithm, 
called the generalized assignment algorithm. Then, we construct a bipartite 



graph G and apply our new algorithm on G. We also improved an 0(n 3 ) 
algorithm for the limited-capacity assignment problem. 
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